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We make use of isospin constraints to study the parametric coupling model and the properties of asymmetric 
nuclear matter. Besides the usual constraints for nuclear matter - effective nucleon mass and the incompress- 
ibility at saturation density - and the neutron star constraints - maximum mass and radius - we have studied the 
properties related with the symmetry energy. These properties have constrained to a small range the parame- 
ters of the model. We have applied our results to study the thermodynamic instabilities in the liquid-gas phase 
transition as well as the neutron star configurations. 



I. INTRODUCTION 

The search for an effective nuclear matter theory has been 
one of the main goals of nuclear physics (see, for instance, the 
pioneering work of Fetter and Walecka [1]). The advances in 
experimental data acquisition has played an essential role on 
the improvement of nuclear models. One of the most required 
parameters to constrain symmetric nuclear matter models, the 
incompressibility [2], has had many contributions since the 
non-linear Boguta and Bodmer model was proposed as an 
alternative to the Walecka model J3l • The derivative coupling 
models proposed by Zimanyi and Moszkowski were also in- 
spired to correct distortions in the incompressibility and effec- 
tive nucleon mass determinations |5|]. 

In this work we make use of the so called parametric cou- 
pling model (PCM) [6] which describes nuclear matter with 
density dependent coupling constants. The PCM, presented 
in the next section has appeared in the literature after a first 
attempt to unify Walecka and ZM models in an unique model 
J3]. The PCM was successful in reproducing the results ob- 
tained in the Walecka [4], ZM [5] and exponential coupling 
models [8] and even more successful in reproducing nuclear 
matter and neutron star main properties using special choices 
for the free parameters that define the density dependence of 
the coupling constants. The PCM has been used to study a 
wide range of problems, specially the nuclear matter com- 
pression modulus |9|], the inclusion of strange meson fields 
[10] and the hadron-quark phase transition jlll,[l2ll . 

Presently, an important issue in nuclear physics is to con- 
strain the nuclear matter equation of state (EOS) from com- 
pact star properties or from laboratory measurements such as 
the ones planned for the next-generation of exotic radioactive 
beam facilities at CSR/China, FAIR/Germany, RIKEN/Japan, 
SPIPvAL2(GANIL)/France or the planned Facility of Rare Iso- 
tope Beams/USA, where the high density behavior of the sym- 
metry energy will be further studied experimentally. While 
saturation properties of symmetric nuclear matter such as the 
saturation density, binding energy and incompressibility are 
quite well settled, properties of asymmetric nuclear matter, 
such as the density dependence of the symmetry energy are 
much less constrained. 
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Many authors have recently tested nuclear models, both 
phenomenological lfl3l - l23ll . variational [IV] or microscopic 
ill 8l l24ll such as the Brueckner-Hartree-Fock formalism, in 
order to impose constraints or verify their compatibility with 
data coming from heavy ion collisions, giant monopole reso- 
nances or isobaric analog states. In particular, recent exper- 
imental constraints from isospin diffusion in heavy-ion col- 
lisions (HIC) give for the slope parameter of the symmetry 
energy the value L = 88±25 MeV and for the isospin incom- 
pressibility coefficient K T = -500 ± 50 MeV CIS. The 
latter is in agreement with the value of K T = —550 ±100 
MeV predicted by the independent measurement of the iso- 
topic dependence of the giant monopole resonance in Sn iso- 
topes HHH] and the value K T = —500 + 1 2 ' 

tained from the study of neutron skins ni 
isoscaling in HIC favors L ~ 65 MeV |2_ 
sion ratios favor L ~ 55 MeV i28tl . 

The main goal of the present work is to constrain even more 
the free parameters in the PCM by studying the symmetry en- 
ergy and its density dependence through its slope and its sec- 
ond derivative with respect to the density. We try to establish 
some correlations between the symmetry energy and its den- 
sity derivatives with quantities with astrophysical interest such 
as the crust-core crossing density, the density at muon and at 
strangeness onset, the maximum star mass, and the threshold 
star mass for direct Urea lf29ll . 

In section|n]we summarize the PCM, in section[nI]the sym- 
metry energy dependence on the density is discussed, the ther- 
modynamic instabilities are discussed in section [TV] and the 
neutron star properties in section [V] Finally we draw some 
conclusions in section [VTl 



n. PARAMETRIC COUPLING MODEL 

We first study the nuclear matter properties within the PCM. 
This model is based on the Walecka model 0]. However, in 
order to correct the values of the effective nucleon mass and 
incompressibility of symmetric nuclear matter at saturation 
density, Taurines et. al. |6|] have introduced parametric cou- 



2 



plings. They are introduced in a similar way as Zimanyi and 
Moszkowski included the derivative couplings in their mod- 
els. Here we briefly describe the model and the reader is ad- 
dressed to Ref. [6] for further details. The PCM Lagrangian 
density reads 

L = - (Mb - g* aB &) - gtsl^^B 

B 

- J2^b[^9*b1^ t ■ P^b + ^4>\[ii^ - m x }ip\ 

B A 

- \p pv ■ p"" + \ m% P P» ■ 

(1) 

where 

9*aB = m * a B9<y ; 9uB = ™*pb9w 5 9p B = m* B g p (2) 

and 

<B = (1+Hrr)~ n ; n = a,p, 7 . (3) 
nlviB 

In the equations above, ipB represents the baryon fields that 
can be summed over the whole baryon octet. The baryon 
fields are coupled to the meson fields a, to and p. The electron 
and muon fields appear as ip\ and must be introduced for the 
description of stellar matter. 

The parametric couplings restore the Walecka model for 
a = f3 = 7 = 0, the ZM1 model for a — 1 and j3 = 7 = 0, 
the ZM3 model for a = /3 = 1 and 7 = 0, and exponential 
coupling model for a = (3 = 7 = 00. 

The 9cr an d 9u coupling constants are chosen to reproduce 
the binding energy E B = e/p — M = -15.75 MeV at the 
saturation density po = 0.16 fm -3 . We fix the isovector cou- 
pling constant, g p to fit the symmetry energy, 04 = 32.5 MeV, 
at the same saturation density. 

We will work with four versions of the PCM: a) varying 
the scalar parameter a with the vector parameter /3 = and 
isovector parameter 7 = (Model-S); b) varying the scalar 
and vector parameters a — f3 while keeping the isovector 
parameter 7 = (Model-SV); c) varying all three parame- 
ters a = f3 = 7 (Model-SVI) and d) varying the scalar and 
isovector parameters a = 7 while keeping the vector parame- 
ter (3 — (Model-SI). We expect to verify the nuclear matter 
main properties dependence on the coupling constants. We 
will also consider the parametrization a) SI2 with a = 0.1, 
(3 = and b) SVI2 a = (3 = 0.2 while, in both cases, we 
vary 7. These last two choices will allow us to explore the 
isovector degree of freedom while choosing reasonable prop- 
erties for the isoscalar channel, namely for symmetric nuclear 
matter. 

In Fig. Q] we represent the meson couplings g a , g u and g p , 
as a function of the parameter a for all models. We plot the 
couplings for models S and SI in one figure and the ones of 
models S V and S VI in the other. Since the values of the g a and 
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FIG. 1: The coupling parameters gi of the models a) S and SI; b) SV 
and SVI, as a function of a. For each couple of models the g a and 
g u parameters are equal and only the g p couplings vary. g a (solid 
line), g^ (dotted line), g p (dashed line for S and SV, dot-dashed line 
for SI. 



g^ couplings are first chosen to reproduce symmetric nuclear 
matter data and only after this choice we establish the value of 
g p , the scalar and vector coupling constants are the same for 
the pairs of models S, SI and SV, SVI. These two couples of 
models differ only in the isovector channel. It is seen that for 
models S and SI the g p parameter has stabilized at a constant 
value for a ~ 0.2. For the models SV and SVI this occurs 
only for a ~ 0.5. The g a and g u parameters stabilize only for 
a ~ 1, however, saturation properties of symmetric nuclear 
matter restrict the accepted values to a < 1. In all four dif- 
ferent versions of the model, a = 00 has a stable result which 
essentially does not differ from a = 2 and the changes from 
a = 1 are very small. The density dependence of the coupling 
constants can be examined from Fig|2j where we have plotted 

2i. = m* = (l + 9 " a \ 
gi m « V nM B J ■ 

The greater the value of n the stronger the density depen- 
dence and the reduction of the effective coupling constant 
value at higher densities. In order to compare the results ob- 
tained with PCM, we have also plotted, in the same figure, the 
density dependence of the relative TW couplings 13011. The 
behavior of r^/r^o, i = cr, u> is similar to the one obtained 
with PCM for a = 1, The same does not occur for T p which 
has a much faster decrease with density. 
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FIG. 2: (Color online) The effective coupling constant to bare 
coupling constant ratio g*/gi for n=0.0 (full), 0.05 (brown, long 
dashed), 0.20 (dark green, medium dashed), 0.40 (pink, short dashed) 
and 1.0 (cyan, dotted line) with i — a, u, p for different choices of a 
in the PCM, as well as TW functionals r a (purple, short dot-dashed), 
r u (light green, long dot-dashed) and r p (light brown, dot-dot-dashed 
line), where r« = I^/To. Notice that TW curves are shown with 
thicker lines. 
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FIG. 3: The dependence of the slope of symmetry energy on the free 
parameter a for model S (dotted line), SV (solid line), SI (dot-dashed 
line) and SVI (dashed line). 



We have included hyperons in the present model for the 
high density EOS in the inner regions of a neutron star. Since 
the discussion of the hyperon coupling constants is not the aim 
of this work, we just fix the hyperon-meson coupling con- 
stants to X = — = — = — = \/273- In order to 

g°N 9ioN 9pN V ' 

study neutron stars we must also include the Baym-Pethick- 
Sutherland (BPS) EOS for densities below the neutron drip 
line Hfll. 
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FIG. 4: (Color online) a) Symmetry energy, b) the symmetry energy 
slope L and c) the incompressibility K 3ym for model SI with a=0, 
0.10, 0.20, 1.00=7; /9=0- legend: a = (full line), 0.1 (dark brown 
long dashed line), 0.2 (green short dashed line), 1.0 (cian dotted line), 
NL3 (purple dot-dashed line) and TW (light brown dot-dot-dashed 
line). 



III. SYMMETRY ENERGY 

In Tables [ BlVI we give some of the properties of 
the models under study: the symmetry energy slope 
L = 3podE sym /dp, the symmetry energy incompressibility 
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FIG. 5: (Color online) K aym and K T as a function of the symmetry energy slope L for all the parametrizations studied. The presently accepted 
limits from isospin diffusion (ID) in heavy-ion collisions 1161 [Tal . giant monopole resonance (GMR) in Sn isotopes [25, 26] and neutron skin 
studies (NS) Hi . 



K sym = 9pod 2 E sym /dp 2 , the symmetry term of the incom- 
pressibility of the nuclear EOS K T — K sym — L(6 — Qo/ K), 
where Qo = 27 p^d 3 (E / A) / dp 3 , the effective nucleon mass 
M* I'M, the onset density of strangeness p s , the onset density 
of muons p„, the threshold density for the direct Urea process 
PDU and the mass of a star with central density equal to Pdu> 
Mdu /Mq, the mass Mm ax /Mq and radius of the stable star 
configuration with maximum mass. 

The slope L is within the experimental constraints 55 MeV 
< L < 113 MeV for all parametrizations considered. In the 
S and SV models the slope has a quite large value, though 
never below 95 MeV. The variation of L as the free parameter 
a increases is shown in Figj3] It is seen that L decreases with 
the increase of a, starting at ~ 110 MeV and stabilizing at 
- 95 MeV for S and SV models, at - 75 MeV for the SI 
model, and at - 60 MeV for SVI. 

A value of K T in the interval [-400, -675] MeV puts very 
strong restrictions. According to some authors ll22l l32i 13311 
however, it is difficult to determine the experimental value of 
K T accurately which may suffer from the same ambiguities 
already encountered in earlier attempts to extract the incom- 
pressibility coefficient of infinite nuclear matter from finite- 
nuclei extrapolations. We, therefore, do not rule out models 
with values of K T larger than -400 MeV. In fact, in (H^l the 
properties of several relativistic and Skyrme models which re- 
produce the ground state properties of stable and unstable nu- 
clei or the properties of nuclear or neutron matter have been 
compared and it was shown that a large number of models 
have K T < -400 MeV. 

For reference, in Fig. [4] we plot for the different 
parametrizations considered the symmetry energy (a), its 



slope L (b) and the incompressibility K sym (c), for model 
SI. We also include the corresponding curves obtained with 
the non-linear Walecka model parametrization NL3 0411 and 
density dependent relativistic hadronic model TW 83011 . Vary- 
ing a and 7 between and 1 , we go from the hard symmetry 
energy behavior of NL3 to the soft TW behavior. Similar con- 
clusions are taken for the slope and the incompressibility: as 
we increase a the behavior of our parametric model changes 
from a NL3-like to a TW-like behavior. 

In order to compare the predictions of the different 
parametrizations under study and the experimental constraints 
on L and K T we plot for all the models K sym and K T versus 
the slope L in Fig. [5] It is seen K sym and K T show a linear 
correlation with L as also discussed in (24]]. We have checked 
that in fact this correlation also exists for the ratio Qo/K 
which enters the definition of K T . We include the constraints 
obtained from isospin diffusion (ID) in heavy-ion collisions 
llrjl [Hil . isotopic behavior of the giant monopole resonance 
(GMR) in Sn isotopes JUIH] and neutron skin studies (NS) 
II19I1 . The models which better satisfy the constraints are: S 
and SI, with 0.05 < a < 0.1 and SVI, with 0.05 < a < 0.2. 
Model SV only presents values with good agreement with L 
and K T for a ~ 0.2. So far we have considered the S, SV, 
SI and SVI versions of the PCM in which two, one or none of 
the parameters has been set to zero, while the remaining ones 
have been set equal. In the next sections we still keep this 
analysis to impose some constraints on the parameter choice, 
using neutron stars observations. We finally vary the param- 
eters independently and choose two sets which give excellent 
values not only for L and K T but also for K and neutron star 
maximum masses. 
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TABLE I: Results for Model-S. The binding energy (15.75 MeV), saturation density (0.16 fm 3 ) and symmetry energy (32.5 MeV) at satu- 
ration are the same for all parametrizations. The quantities given are the symmetry energy slope L, the symmetry energy incompressibility 
K aym , the symmetry term of the incompressibility of the nuclear EOS K T , the effective nucleon mass at saturation density M* / M, the onset 
density of strangeness p s , the onset density of muons p M , the threshold density for the direct Urea process pou and the mass of a star with 
central density equal to pou, Mou /Mq, the mass Mmax /Mq and radius of the stable star configuration with maximum mass. 
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TABLE II: The same of Table|I]for Model-SI 
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FIG. 6: (Color online) Spinodal curves for S, a = (short dashed), 
S and SI, a = 1 (full and brown dot-dashed), SV and SVI a = 1 
(long dashed and cyan, dotted lines), as indicated. The red, dot- 
dashed and dark blue dotted lines represent the /3-equilibrium EOS 
for a — /3 = 7 = 1 (red) and a = /3 = 7 = (blue line), and the 
circle, triangle and square stress the crossing point of the EOS with 
the spinodal for their respective a values. The inset magnifies the 
crossing of the /3-equilibrium EOS with the spinodal. 



IV. THERMODYNAMICAL INSTABILITIES 

The system is unstable against phase separation keeping 
volume and temperature constant, if the free energy curvature 
matrix J 7 ^ 

_ / d 2 T \ 
13 \dpidpJ T 

is negative 13511 . Stability implies that the free energy den- 
sity J 7 is a convex function of the densities p p and p n . In the 
present section, we define the free energy T with no leptons 
or hyperons. 

Since we are interested in studying the spinodal curve of 
the two-fluid nuclear system, it is enough to evaluate the zero 
of the smallest eigenvalue A_ of J 7 ^ , where 

A_ = i (Tr(7) ± v/Tr(^) 2 - 4Det(^)) , 
and where the relevant eigenvectors are 

Spj~ _ A_ — Fa 
8 pi" Fij 

wifhi,j =p,n Sill]. 

The thermodynamic spinodal, defined by the A_ = con- 
dition determines the instability regions of the system and is 
plotted in Fig. [6] for the models under discussion. We point 
out that, for symmetric matter, all the models have the same 
density, binding energy and symmetry energy at saturation 
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TABLE III: The same of TableUfor Model-S V 





L 


sym 


K 


K T 




A/f * / A/f 
1V1 / ivl 


Pt 


Pa 




Pdu 


A/f 

Mdu 


A/f 

Mm ax 


K 




(MeV) (MeV) (MeV) (MeV) 


^iviev ) 




(rm ) 


(1m ) 


(im ) 


(1m ) 


(A/f \ 


( A/f \ 
(Mq) 


(™) 


0.00 


111.42 


105.67 


566 


-970 


-2068.27 


0.537 


0.1025 


0.244 


0.109 


0.194 


0.77 


2.58 


12.61 


0.05 


108.29 


83.26 


458 


-848 


-1190.65 


0.554 


0.0938 


0.254 


0.108 


0.197 


0.77 


2.42 


11.96 


0.10 


105.84 


63.78 


360 


-751 


-611.36 


0.570 


0.0861 


0.263 


0.108 


0.199 


0.73 


2.26 


11.31 


0.20 


102.40 


36.71 


276 


-560 


47.68 


0.600 


0.0793 


0.281 


0.108 


0.204 


0.69 


1.96 


10.06 


0.40 


99.00 


12.42 


195 


-336 


483.72 


0.649 


0.0727 


0.310 


0.108 


0.212 


0.63 


1.39 


8.70 


1.00 


96.66 


-1.99 


159 


-165 


685.86 


0.710 


0.0699 


0.342 


0.109 


0.221 


0.62 


0.83 


9.00 



TABLE IV: The same of TableUfor Model-SVI 
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but different incompressibilities and effective nucleon masses. 
The model with the largest incompressibility (S, a = 0) has 
the largest spinodal and the ones with the smallest incom- 
pressibility (SV and SVI, a = 1) have the smallest spinodal. 
The point of the spinodal with p p = p n corresponds to the 
minimum of the pressure of symmetric nuclear matter, when 
the incompressibility is zero. Since we impose that the pres- 
sure is zero for all models at the same density, the saturation 
density is 0.16 fm~ 3 , the incompressibility of the models with 
a larger curvature will become zero at a larger density. We no- 
tice that the effect of the parameter 7 is very small and only 
seen for very asymmetric matter, when it slightly reduces the 
spinodal. 

In Fig. [6] we also plot the /3-equilibrium EOS for a = 
and for model SVI with a — 1. The crossing of the (3- 
equilibrium EOS with the thermodynamical spinodal instabil- 
ity line gives a prediction for the transition density approxi- 
mately 15% larger than the value obtained from a Thomas- 
Fermi calculation of the pasta phase 13711 . Therefore, we de- 
termine the crust-core transition density from the crossing be- 
tween the EOS and the spinodal and expect that our estimation 
of the transition density will define an upper bound of the cor- 
rect transition density 113 811 between the crust and the core of 
a neutron star. It is seen that quite different predictions are 
obtained within PCM for the transition density both because 
the spinodals have a different behavior for neutron rich matter 
and the EOS for /3-equilibrium stellar matter also differ when 
varying a from to 1 . Even if the EOS do not differ consid- 
erably for different parameters, the values of the density and 
isospin asymmetries at the crossing of the EOS with the spin- 
odal are found to be quite different, though always with low 
proton densities. In Fig. [6] we only plot the EOS for a = 
and a = 1 for SVI. We have obtained for the transition den- 
sities: 0.102 fm" 3 (a = 0), p t ~ 0.095 fm" 3 (models S and 



SI with a = 1) and p t ~ 0.070 fm" 3 (models SV and SVI 
with a = 1). The values of p t given in Tables HKlIVI show 
a decrease of p t with the decrease of the slope L, contrary to 
the results of 112411 . However, we have to take these values with 
care because we are not only varying the slope L but also the 
incompressibility K while imposing the same saturation den- 
sity for all models. As a consequence, the extension of the 
spinodal region grows with K and it is this correlation that 
is reflected in the data of Tables Ul lIVI In the next section we 
calculate the transition density for a selected set of parameters 
and then we will be able to obtain a correlation between L and 
p t similar to the one obtained in 12411 . 



V. NEUTRON STAR CONSTRAINTS 

Besides heavy ion experiments, we can make use of astro- 
physical observations of neutron stars to stablish some con- 
straints on the PCM and its free parameters. We therefore ob- 
tain the EOS of /3-equilibrium matter, including the presence 
of hyperons, and integrate the Tolman-Oppenheimer- Volkoff 
equation for a spherical compact star in hydrostatic equilib- 
rium 091 kfJl in order to obtain the family of stars which cor- 
respond to each parametrization considered. 

The different EOS are shown in Fig|7^) for all four versions 
of the model, considering the limiting values a = 0.05 and 
a = 1.00. The results for values of 0.05 < a < 1.00 are 
intermediate and lie between the two curves. The kink on each 
curve corresponds to the onset of strangeness. For a = 0.05 
the incompressibility is 3 10 MeV for S and SI models and 458 
MeV for SV and SVI. The last EOS is the hardest one and 
therefore, the onset of hyperons occurs at smaller densities 
for this model. The hyperons make the EOS softer but never 
softer than any of the other versions (with higher values of a) 
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FIG. 7: (Color online) a) Equation of state and b) it correspond- 
ing Mass - Radius curves for the families of stars obtained with 
the parametrizations of Tables IHIVI As indicated on the figure, for 
a=0.05 the coincident S and SI curves are shown by the green, short 
dashed line, whereas the red, dot-dot-dashed line stands for the SV, 
SVI single curve. For a=1.0, the S, SI, SV, SVI curves are shown by 
full, brown dot-dashed, long dashed and cyan dotted lines, respec- 
tively. 



of the model. 

The corresponding neutron star families (mass-radius) are 
shown in Fig|7p). As expected, one can see that all mod- 
els with low values of a have maximum neutron star masses 
above 2M Q , while the results drastically change when val- 
ues of a approach the unit, reducing the maximum neutron 
star mass obtained with the S and SI models to values close to 
1 .5Mq. The results get even smaller when S V and SVI model 
are considered, reducing the maximum neutron star mass to 
values below 1 Mq . 

One can stablish the minimum value required for the max- 
imum neutron star mass from observations. From Tables HlUVl 
we see that, for the S and SI models, which present almost the 
same neutron star masses, a < 0.1 allow neutron star masses 
greater the 1 .8Af Q , while the models S V and SVI have a wider 
range of values (a < 0.3) that yield the same star mass. 



According to 14111 . observations on cooling predict that 
direct Urea should occur only in stars with a mass M > 
1.35M . We have determined the star with maximum mass 
that does not allow the direct Urea process {Mdu) f° r the dif- 
ferent parametrizations of the model proposed. By pr>u we 
denote the baryonic density that defines the onset of the Urea 
process. The results are summarized in Tables HlUVl No model 
satisfies the condition Mdu > 1.35Mq. All four versions 
of the model predict the Urea process with densities 0.194 
fm~ 3 < pdu < 0.293 fm~ 3 which correspond to neutron 
stars masses around 0.57 Mq < Mdu < 0.77 A/ Q . This re- 
sult is expected as PCM presents the same main features as 
the models discussed in 14111 which are not in agreement with 
Urea expectations. The threshold density for the direct Urea 
process is plotted versus the slope L in Fig. [8] It is seen that 
it strongly depends on the slope of the symmetry energy: a 
smaller slope corresponds to a larger density of the Urea on- 
set. 

We expect to be able to infer isospin asymmetry properties 
or stablish a few constraints on the parametrizations from neu- 
tron star observations. When comparing S and SI models or 
SV and SVI models, which present the same global behavior, 
we have identified a small difference in the neutron star radii 
that becomes larger as a increases. For the S and SI mod- 
els the difference between the radii reaches Ai?=380 m, for 
a = 1.0 while it reaches Ai?=330 m for SV and SVI model. 
This difference would correspond, for example, to a differ- 
ence of Avf ~ 0.1 kHz (~ 4%) in the fundamental mode of 
a star i42l 14311 and is larger if we consider stars with masses 
below the maximum allowed mass. 

In Fig |9] several properties of the neutron star structure are 
plotted as a function of the symmetry energy slope for the 
models S, S V, SI and SVI. The figure shows the baryonic den- 
sities at the onset of hyperons for each model. Hyperons are 
very important because their appearance results in the soften- 
ing of the EOS and the lowering of the neutron star maximum 
mass. In the same figure, it is also represented the baryonic 
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FIG. 8: Threshold density for the direct Urea process to occur as 
versus the slope L. The cyan dotted, brown dot-dashed, dashed and 
full lines stand for the SVI, SI, SV and S models, respectively. 
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TABLE V: Results for a = 0.1 and f3 = 

7 M*/M if L K 3ym K T Qo p t p s p M p DU M DU M M ax R 

(MeV) (MeV) (MeV) (MeV) (MeV) (fm~ 3 ) (fm~ 3 ) (fm~ 3 ) (fm~ 3 ) (M ) (M ) (Km) 

0.450 0.749 224 60 -123.5 -585 -378.93 0.0894 0.305 0.107 0.253 0.68 1.811 9.512 

0.240 0.749 224 70 -87.7 -587 -253.76 0.0874 0.305 0.107 0.240 0.67 1.806 9.510 

0.120 0.749 224 80 -46.3 -583 -158.76 0.0865 0.305 0.107 0.231 0.66 1.803 9.511 



TABLE VI: Results f or a = = 0.2 
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FIG. 9: Density for the strangeness (full), the Urea process (dotted) and the muon (dashed line) onset for each version of the model. 



density at the muon onset and the density at which the proton 
fraction reaches the value which allows the direct Urea pro- 
cess to occur in neutron stars (for details see Ref. lr29[ EUd . 
While the muon onset is not very sensitive to the value of L, 
the hyperon onset may vary 40-50%, and the densities for the 
onset of Urea, 20-50% between the limits of the plotted values 
of L which are all within the experimental constrained values. 

In Fig{lO]we analyse the dependence of the maximum mass 
neutron star radius on a and the asymmetry parameters L, 
K sym and K T . It is shown that there is a correlation which 
may help to impose restrictions on the parameters: smaller 
values of L, K sym and \K T \ give larger radius for the maxi- 
mum mass configuration. 

We propose to use the new parametrizations SI2 and SVI2 
because models S, SI, SV and SVI with reasonable incom- 



pressibilities predict too small maximum star masses. With 
a fixed a and/or a — (3 parameter, and changing the 7 pa- 
rameter we have changed the symmetry energy slope. Tables 
IVl and [VI] show the properties of the EOS and neutron stars 
for these new parametrizations. The maximum star masses 
have improved and are of the order 1 .8 M & while the proper- 
ties of symmetric nuclear matter are all reasonable. However 
the density for Urea onset is still too small in all models and 
the radius of the maximum mass configuration is also quite 
small. We have tried to vary all the parameters independently 
of each other but we were not able to find a parametrization 
which predicts larger densities for the onset of the Urea pro- 
cess than the ones already obtained. This seems to be a weak 
point of these models. We made also an estimation of the 
crust-core transition densities by calculating the crossing of 
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the /3-equilibrium EOS with the thermodynamical spinodal as 
explained in section IV. The results are given in Tables [Viand 
I VII and vary between 0.086 and 0.089 fm~ 3 showing a small 
decrease when the slope of the symmetry energy increases. 

We next compare our results with two RMF models: NL3 
with constant couplings and TW with density dependent cou- 
pling parameters (Table IVTIl ). We have obtained for the mod- 
els SI2 and SVI2 better values for the slope of the symme- 
try energy, well within the experimental range, than NL3 and 
TW give. For these two models the slope of the symmetry 
energy corresponds, respectively, to the upper and lower limit 
of L proposed from the experimental values. We have also 
obtained results for K T closer to the experimental ones than 
those found within NL3 and TW, keeping good results for 
the compression modulus, effective nucleon mass and neutron 
star global properties. PCM introduces a density dependence 
on the coupling parameters weaker than the one of TW (see 
Fig. I3, mainly for the isospin channel. The parametrization 
of the coupling parameters within PCM does not allow for 
such a strong density dependence for the coupling parame- 
ters as TW shows. For the crust-core transition density the 
models SI2 and SVI2 preview a density similar to the one ob- 
tained with TW. Concerning the correlation between_p t and 
L, models SI2 and SVI2 behave like those in Ref.||2J, de- 
creasing the value of p t as L increases. The density for the 
strangeness onset is closely related with the incompressibility 
of the EOS: SI2 is very soft and predicts the largest density for 
the strangeness onset, SVI2 gives a result very similar to TW, 
~ 40% larger than NL3. Other properties of SI2 and S VI2 are 
also closer to TW than to NL3 except for the the onset density 
for the Urea process and the mass of the compact star with this 
density at the centre: it is seen that for other similar saturation 
properties, the value of K T for TW seems to be the only great 
difference. 



VI. CONCLUSIONS 

We have employed the parametric coupling model, devel- 
oped by Taurines et. al. Ja], to describe nuclear matter prop- 
erties, in special, those related to the isospin channel which 
were still not investigated in the PCM model. Some of the pa- 
rameters of the model were fixed by the binding energy (15.75 
MeV) and saturation density (0.16 fm~ 3 ) for symmetric mat- 
ter and by the symmetry energy (32.5 MeV). We have tried to 
restrict the remaining parameters by reproducing experimen- 
tal results obtained with heavy ion collisions at intermediate 
energies and neutron star properties. 

Considering a wide range of experimental data for the slope 
of the symmetry energy L, 55 MeV < L < 115 MeV and 
for the symmetry term of the incompressibility of the nuclear 
EOS K T , -375 MeV< K r < -650 MeV, we have ruled out 
some values of the free parameters. The models which best 
satisfy the constraints are: S and SI, with 0.05 < a < 0.1 and 
SVI, with 0.05 < a < 0.2. Model SV only presents values 
with good agreement with the L and K T experimental values 
for 0.2 < a < 0.3. 

Some properties of neutron star matter and neutron star 



structure have been described with the PCM for the range of 
the parameters restricted by the values of L and K T , namely 
the maximum mass star configuration, and the densities of the 
onset of muons, strangeness and the direct Urea process. It 
was shown that all properties were reasonably well described 
except the prediction of the direct Urea process at a too low 
density. It was also shown that the density at the onset of 
muons is almost independent of the parameters of the model 
but the strangeness onset is sensitive to the symmetry energy 
slope: a larger slope corresponds to a smaller onset density. 
The same is true for the direct Urea process, although with 
a slightly weaker dependence. It was also shown that the ra- 
dius and mass of the maximum mass stable star configuration 
was correlated with the L value: a larger L corresponds to a 
larger radius and a larger mass. A similar correlation exists 
with the symmetry term of the incompressibility of nuclear 
matter. However from model to model this correlation has a 
different slope and a given value of L and K T does not define 
in a unique way the neutron star properties. 

Even though the proposed range of validity of the parame- 
ters of the PCM are in agreement with the experimental range 
described above for heavy ions, it is of common belief that 
the slope of symmetry energy should be close to the lower 
limit of the interval 55 MeV< L < 115 MeV JH]. In or- 
der to investigate other possibilities of the model, we have 
proposed new parametrizations chosen to reproduce the sym- 
metry energy slope L=60, 70 and 80 MeV. These parametriza- 
tions have shown excellent results also for K T , K, M* and the 
neutron star maximum mass but still predict quite low densi- 
ties for the onset of the Urea process. The results shown in 
tableslVlandlVIIdescribe some of these results. We have tried 
to find a set of parameters which could predict a larger density 
for the onset of the direct Urea process without success. The 
PC M p resent the same main features of the models discussed 
in 14111 which are not in agreement with Urea expectations. 

We have constrained the free parameters of the PCM using 
the properties of asymmetric nuclear matter, but it is important 
to mention that the EoS which show the best results (tables V 
and VI), for these properties are in perfect agreement with 
those imposed by the analysis of symmetric nuclear matter 
properties such as those described by B44I1 where constraints, 
obtained from the analysis of nuclear matter flow in heavy ion 
collisions, were proposed for the high density EoS of sym- 
metric matter. This problem has been recently examined in 
Ref. j45ll . It is also important to mention that the constraints 
imposed by asymmetric nuclear matter properties are stronger 
than the ones obtained from symmetric nuclear matter (for de- 
tails see Ref.iH). 

We have compared our results with two RMF models: NL3 
with constant couplings and TW with density dependent cou- 
pling parameters (Table fVTTb - It was shown that the models 
SI2 and SVI2 predict results closer to TW except for the on- 
set density of the Urea process. These could be related to 
the larger K T value that TW has since all the other saturation 
properties are similar. It was also shown that while the density 
dependence of the g* and <j>* coupling parameters in PCM is 
similar to the one of TW, the density dependence of g* is very 
different. 
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TABLE VII: Results for TW and NL3 

Model B/A p K M* /M £ sym L K sym K T p t p s p M p DU M DU M M ax R 
(MeV) (fm^ 3 ) (MeV) (MeV) (MeV) (MeV) (MeV) (fm~ 3 ) (fm" 3 ) (fm" 3 ) (fm" 3 ) (M ) (M ) (Km) 

TW 16.3 0.153 240.1 0.56 32.0 55.3 -124.7 -332.1 0.085 0.287 0.115 0.315 1.46 2.012 11.60 
NL3 16.3 0.148 270 0.60 37.4 118.5 100.9 -698.4 0.065 0.217 0.111 0.205 1.00 1.707 12.65 

a L (MeV) 
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FIG. 10: (Color online) Radius of star with maximum mass versus a, L, K sy m and K T , for the four different versions of the PCM: S (full), SI 
(brown dot-dashed), SV (dashed) and SVI (cyan dotted line). 



The relation between the symmetry energy and some prop- 
erties of neutron stars was also discussed. The different mod- 
els considered in the present study with one or more density 
dependent couplings, showed quite large differences among 
them, even for the same values of L and K T . We conclude 
that the properties of nuclear matter for both symmetric and 
asymmetric matter at saturation are not enough to define the 
structure of compact stars and experimental data at two, three 
times saturation density are required. 
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